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Abstract 

This paper presents a statistical framework for assessing wireless systems performance using hierarchical data 
mining techniques. We consider WCDMA (wideband code division multiple access) systems with two-branch 
STTD (space time transmit diversity) and 1/2 rate convolutional coding (forward error correction codes). Monte 
Carlo simulation estimates the bit error probability (BEP) of the system across a wide range of signal-to-noise 
ratios (SNRs). A performance database of simulation runs is collected over a targeted space of system configu- 
rations. This database is then mined to obtain regions of the configuration space that exhibit acceptable average 
performance. The shape of the mined regions illustrates the joint influence of configuration parameters on system 
performance. The role of data mining in this application is to provide explainable and statistically valid design 
conclusions. The research issue is to define statistically meaningful aggregation of data in a manner that permits 
efficient and effective data mining algorithms. We achieve a good compromise between these goals and help 
establish the applicability of data mining for characterizing wireless systems performance. 



1 Introduction 



Data mining is becoming increasingly relevant in simulation methodology and computational science [ RG01 ]. It 
entails the 'non-trivial process of identifying valid, novel, potentially useful, and ultimately understandable patterns 
in data' [ FPSS96 ]. Data mining can be used in both predictive (e.g., quantitative assessment of factors on some 
performance metric) and descriptive (e.g., summarization and system characterization) settings. Our goal in this 
paper is to demonstrate a hierarchical data mining framework applied to the problem of characterizing wireless 
system performance. 

This work is done in the context of the S 4 W problem solving environment [ VHW02 ] — 'Site-Specific System 
Simulator for Wireless System Design'. S 4 W provides site-specific (deterministic) electromagnetic propagation 
models as well as stochastic wireless system models for predicting the performance of wireless systems in specific 
environments, such as office buildings. S 4 W is also designed to support the inclusion of new models into the system, 
visualization of results produced by the models, integration of optimization loops around the models, validation 
of models by comparison with field measurements, and management of the results produced by a large series of 
experiments. In this paper, we study the effect of configuration parameters on the bit error probability (BEP) of a 
system simulated in S 4 W. 

The approach we take is to accumulate a performance database of simulation runs that sweep over a targeted 
space of system configurations. This database is then mined to obtain regions of the configuration space that exhibit 
acceptable average performance. Exploiting prior knowledge about the underlying simulation, organizing the com- 
putational steps in data mining, and interpreting the results at every stage, are important research issues. In addition, 
we bring out the often prevailing tension between making statistically meaningful conclusions and the assumptions 
required for efficient and effective data mining algorithms. This interplay leads to a novel set of problems that we 
address in the context of the wireless systems performance domain. 

Data mining algorithms work in a variety of ways but, for the purposes of this paper, it is helpful to think of them 
as performing systematic aggregation and redescription of data into higher-level objects. Our work can be viewed 
as employing three such layers of aggregation: points, buckets, and regions. Points (configurations) are records in 
the performance database. These records contain configuration parameters as well as unbiased estimates of bit error 
probabilities that we use as performance metrics. Buckets represent averages of points. We use buckets to reduce 
data dimensionality to two, which is the most convenient number of dimensions for visualization. Finally, buckets 
are aggregated into 2D regions of constrained shape. We find regions of buckets where we are most confident that the 
configurations exhibit acceptable average performance. The shapes of these regions illustrate the nature of the joint 
influence of the two selected configuration parameters on the configuration performance. Specific region attributes, 
such as region width, provide estimates for the thresholds of sensitivity of configurations to variations in parameter 
values. 



1.1 Reader's Guide 



Our major contribution is the development of a statistical framework for assessing wireless system performance using 
data mining techniques. The following section outlines wireless systems performance simulation methodology and 
develops a statistical framework for spatial aggregation of simulation results. Section || demonstrates a substantial 
subset of this framework in the context of a performance study of WCDMA (wideband code division multiple 
access [HTOC]) systems that employ two-branch STTD (space-time transmit diversity [ Ala98| ]) techniques and 1/2 
rate convolutional coding (forward error correction codes [ HTOC ]). We study the effect of power imbalance between 
the branches on the BEP of the system across a wide range of average signal-to-noise ratios (SNRs). Section |] 
extends the statistical framework to support computation of optimized regions of the bucket space. Such regions 
are computed by a well-known data mining algorithm [ FMMTO T] , |YFM97 ]. Section || applies these concepts to the 
example in Section |3[ Section |6| summarizes present findings and outlines directions for future research. 
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Figure 1: Typical ID slices of the configuration space. The plots show simulated BERs (bit error rates) of wireless 



systems for five common benchmark channels [ UMT97 ] across a typical range of average SNRs. 



2 The Statistics of Aggregation and the Aggregation of Statistics 



Temporal variations in wireless channels have been extensively studied in the literature [ZT02]. The present work 
uses a Monte Carlo simulation of WCDMA wireless systems to study the effect of these variations. The simula- 
tion traces a number of frames of random information bits through the encoding filters, the channel (a Rayleigh 



fading linear filter [Has93]), and the decoding filters. The inputs are hardware parameters, average SNR, channel 
impulse response, and the number of frames to simulate. The output is the bit error rate — the ratio of the number 
of information bits decoded in error to the total number of information bits simulated. Simulations of this kind 
statistically model channel variations due to changes in the environment and device movement across a small geo- 
graphical area {small-scale fading [ |Has93 ]). We refer to this kind of channel variation as temporal variation because 



a system is simulated over a period of time. Further, we say that a given list of inputs to the WCDMA simulation is 
a configuration or a point in the configuration space. 

Spatial variations are due to changes in system configurations. We use this term to describe two quite different 
phenomena: changes in the average SNR and channel impulse response due to large-scale fading pas93 ] and 



variations of hardware parameters. A typical approach to the analysis of spatial variations is to run several temporal 
variation simulations (i.e., compute bit error rates — BERs — at several points within a given area of interest) and 
plot ID or 2D slices of the configuration space, as shown in Figure [T[ In this paper, we augment this approach with 
statistically meaningful aggregation of performance estimates across several points. The result of this aggregation 
is a space of buckets, each bucket representing the aggregation of a number of points. Moving up one level of 
aggregation in this manner allows us to bring data mining algorithms to operate at the level of buckets. The space 
of buckets mined by the data mining algorithm is then visualized using color maps. The color of each bucket is the 
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c, C, R entities (points, buckets, regions) 

x,b,B random variables 

E[x], E[b], E[B] true means of random variables x,b, B 

a 2 , E 2 true variances of random variables b, B 

x,b,B estimates of means E[x], E[b], E[B] of random variables x,b, B 

a 2 , E 2 estimates of variances a 2 , E 2 of random variables b, B 

P(E) probability of event E, where E is a boolean condition 

F N -i (T) P(X < T) for X having the Student t distribution with N — 1 degrees of freedom 

i x k}k=l set {xi,X2, ••• ,x n } 

{{Zfcj}j = l}fc =1 Set {xn,Xi2, ■ ■ ■ , X lni , X 2 l, X 2 2, ■■ ■ ,X 2rl2 ,... ,X n l,X n2 ,... , X nrln } 



Table 1: Summary of mathematical notation. Lower case letters are used for points and upper case letters are used 
for buckets and regions. Additional conventions are introduced in Table ^ 



confidence that the points (configurations) that map to this bucket exhibit acceptable average performance. 



2.1 The First Level of Aggregation: Points 

Table [I] summarizes some of the syntactic conventions used in this paper. Mathematically, we can think of the 



WCDMA simulation as estimating the mean E[xk] of a random variable xj~ with some (unknown) distribution [jfBS92|] 
(Xk is one when the information bit is decoded in error or zero when it is decoded correctly). Each BER x^j, 
1 < j < nk, output by the simulation is an unbiased estimate of the BEP E[xk] of the simulated configuration Cfc. 
Instead of building a detailed stochastic model of the simulation (analytically, from the distribution of x^), we choose 
to work with the simpler distribution of the BER Xkj, referred to henceforth as just bk- Thus, each sample from the 
distribution of bk is realized by simulating a number of frames and obtaining an estimate of -Efx^]. The distribution 
of bk is approximately Gaussian due to the Central Limit Theorem. Technically, we assume that the number of 
frames per estimate Xkj is 'large enough' so that the Lindeberg condition is satisfied, that the variance of x^ is 
finite, and that {xkjVj=i are i.i.d. We say that E[bk] = E[E[xk]] is the expected BEP of configuration cu under 
Rayleigh fading. 



2.2 The Second Level of Aggregation: Buckets 

Let us now aggregate several points (i.e., random variables) into one bucket. The purpose of this aggregation is to 
reduce data dimensionality to a size that is easy to visualize, usually one or two dimensions. The basic idea is to 
linearly average all points that map to the same bucket but we must do so carefully, in order to preserve a meaningful 
statistical interpretation. Let {&fc}£ =1 be Gaussian random variables with means {-E[bfc]}£ =1 and variances {cr^}^ =1 . 
As in the previous paragraph, let each such variable bk be the estimated BEP of some configuration Ck, 1 < k < n. 
For bucket C, define a bucket (mixture) random variable B as the convex combination 

n 

B = y^Pkbk, 
k=l 

where the pk > and J2k=iPk = !• It is convenient to make {pk}k=i tne probabilities of occurrence of the 
configurations {cfc}£ =1 in the dataset being analyzed. This setup underlines the dependence of the outputs on the 
distribution of the inputs and frees the user from having to provide values for the constants {pkYk=v It i s wei l known 
that, as long as {bk}k=i are mutually independent and Gaussian with means {^[frfc]}£ =1 and variances {cr|}^ =1 , B 
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is Gaussian with mean E[B] = Y^k=xPkE\bk] an(1 variance £ 2 = Y^k=\Pk a \ QCB02| ]. The expected value E[B] 
of the random variable B can be viewed as the expected BEP of bucket C = {c\, C2, • • • , c n } in a Rayleigh fading 
environment, conditional on the (discrete) distribution of the configurations in C. 

The values {pk}k=i are what the statisticians call prior probabilities. For most purposes of this paper, we simply 
estimate {pk}k=i from available data. These values are explicitly or implicitly constructed during experiment design 
and we assume that they remain constant during experiment analysis. However, one can collect additional data as 
long as doing so does not change {pk}k=i- Prior probabilities can come from a number of sources: channel sounding 
measurements, propagation simulations, hardware and budget constraints, or even educated guesses by wireless 
system designers. The rest of the paper silently assumes that the values {pk}k=i have been established beforehand. 
It is important to remember that even though the prior probabilities are for the most part transparent to the analysis 
presented here, they nonetheless always exist and all conclusions of data analysis are made conditional on the prior 
probabilities. 

This discussion of {pk}k=i can be interpreted as a deferral of the exact definition of B until experiment setup, 
or as parameterization of the analysis procedure. A natural question is whether or not this level of parameterization 
is sufficient. It is sufficient for the purposes of this paper but, strictly speaking, the interrelations between {&fc}£ =1 
should also be defined during experiment setup. Mutual independence of {6fc}^ =1 is a simplifying assumption and 
it might be desirable to model interactions between {6^}^ =1 in practice. This implies adding covariance terms to E 2 
and re-thinking the distribution of B. Such analysis is necessarily specific to a particular experiment. For the sake 
of simplicity, the rest of this paper assumes mutual independence of variables in a given bucket. 



2.3 Confidence Estimation 

Point and bucket estimates of the expected BEP are meaningful performance metrics for wireless systems. Let us 
also estimate our confidence in these estimates. Confidence analysis enables wireless system designers to make 
more practical claims than point estimates alone. A statement of the form 'this configuration will exhibit acceptable 
performance in 95% of the cases' is often preferable to a statement of the form 'the expected BEP of this config- 
uration is approximately 5 x 10~ 4 '. More precisely, we say that configuration exhibits acceptable performance 
when the expected BEP E[bk] of configuration is below some fixed threshold T. This statement is conditional 
on the temporal simulation assumptions, i.e., Rayleigh fading. Standard values for T are 10" 3 for voice quality 
systems and 10~ 6 for data quality systems. Likewise, we say that bucket C (a subspace of configurations) exhibits 
acceptable average performance when the expected BEP E[B] of bucket C is below some fixed threshold T. This 
statement is conditional on both the temporal simulation assumptions and the distribution of configurations {c^}^ =1 
in the bucket (the prior probabilities). 

The confidence that configuration (resp. bucket C) exhibits acceptable (average) performance is P(E[bk] < 
T) (resp. P(E[B] < T)). Since b^ and B are Gaussian, these probabilities can be estimated as 

P(E lh] < T) « (1^) , P(E [B] < T) - Fjv-i (0J , 

where Fjy-i (") is the CDF of the Student t distribution with N — 1 degrees of freedom and and N are the sample 
sizes for configuration and bucket C, respectively. For configuration c&, 

^ = - E %' % = j—T) E(% - ^) 2 ' 

where b\ and o\ are the estimates of the expected BEP and the BEP variance at point Cfc, > 2 is sample size, 
and {xkj}™=i are sample values. For bucket C, we substitute point estimates into E[B] = Y^k=iPkE\bk] and 
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£ 2 = £Li^ 2 to obtain 
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where and E are the estimates of the expected BEP and the BEP variance at bucket C, and {pk}^ = i are the prior 
probabilities estimated from the dataset as p^ = n^j Ya=i n » - Observe that 



n ^ n n k 

b = Y^pkh = v » - Yl Yl *u 

k =i ^ k=i nk k=i j=i 

is exactly the sample mean of all observations in the bucket, but E 2 is not the variance of this sample. This is the 
case because {{xkj}^t\}k=i we not samples from the mixture distribution of B — they are samples from the 
constituent distributions of {frfc}£ =1 - 



3 Extended Example 

Let us now apply the techniques developed so far to analyze the performance of a space of configurations. The wire- 
less systems under consideration employ WCDMA technology with two-branch STTD and 1/2 rate convolutional 
coding. We require that the transmitter has two antennas (branches) separated by a distance large enough for their 
signals to be uncorrected, but small enough for the mean path losses and impulse responses of their channels to be 
approximately equal at receiver locations of interest. We assume Rayleigh flat fading channels, which is reasonable 
for indoor applications in the ISM and UNII carrier frequency bands (2.4 and 5.2 GHz, respectively). The goal is to 
study the effect of power imbalance between the branches on the BEP of the configurations across a wide range of 
average SNRs. 

This section presents a number of plots that summarize simulated BERs. We also outline the process of statis- 
tically significant sampling of the configuration space. The next section develops a data mining methodology that 
solves a practically important problem: given a dataset similar to the one presented next, find a region of the config- 
uration space where we can confidently claim that configurations will exhibit acceptable (average) performance. 

Let us begin with an initial sample of the configuration space, as shown in Figure || (top). This figure shows the 
simulated BER as a 2D function /(Si, S 2 ) of the average branch bit energy-to-noise ratios (SNRs) Si and S 2 , in dB. 
The parallel simulation ran for three days on 120 machines (AMD Athlon 1.0 GHz) at a speed of approximately 2.5 
points per machine per day. 10000 frames, or 800000 information bits, were simulated for each of the 820 points 
S 2 = 3, 4, . . . , 42; 5i = 3, 4, . . . , S 2 . Since f(S u S 2 ) is symmetric JDDJS02| ], we show M = 1600 points {c k }jf =1 
for a full cross-product of Si and S 2 . 

Wireless system designers are more accustomed to ID slices of the configuration space, e.g., the ones shown in 
Figure |3[ Define the branch power imbalance factor 

a =10-0-1151-521, 

where S\ and S 2 are the average SNRs of the branches, in dB. (This definition applies as long as the mean path 
losses of the branches are equal.) By definition, < a < 1, where zero corresponds to a total malfunction of one of 
the branches and one corresponds to a perfect balance of branch powers. The graphs in Figure || were obtained by 
fixing a and varying the effective SNR 

S = 101og 10 ((10 - 151 + 10 152 )/2) , 

in dB (top), and fixing the effective SNR and varying a (bottom). (Note that fixing the effective SNR is equivalent 
to fixing total transmitter power.) The sample of configurations came from the dataset shown in Figure || (bottom), 
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Figure 2: (top) Estimates of the BEPs for a space of configurations {ck}jf =1 (M = 1600 points at 10000 frames 
per point). The X and Y axes are the average SNRs of the branches (in dB). The Z axis is the (base ten) logarithm 
of the simulated BER. These estimates are not statistically significant, (bottom) Statistically significant estimates 
{bk}^L 1 of the expected BEPs {E[b/,}}ff =1 for the same space of configurations {ck}%L v For the most part, we are 
90% confident that the estimated expected BEP lies within 10% of its true value. See text for exceptions. 




Figure 3: ID slices of the configuration space {cfc}^£ 1 with fixed branch power imbalance factor a = 10 01 l Sl 52 
and varying effective SNR S = 101og 10 ((lO 0151 + 10 1S2 )/2) (top), and fixed effective SNR S and varying 
branch power imbalance factor a (bottom). These slices were computed from the surface fit onto the data in Figure || 
(bottom). The entire fitted surface is shown in Figure |J 
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Figure 4: A surface fitted onto the statistically significant results in Figure || (bottom). We used a local linear least 
squares regression with a 5% neighborhood and tricubic weighting. This procedure was chosen because it can 
approximate the relatively steep edge of the tolerance region. See [HTF01] for details. 



described in detail later. However, this sample does not contain the exact points for typical slices, so we used a fitted 
surface — Figure || — to approximate the BERs for the slices in Figure |3| We choose to work with the axes 5*1, S2 in 
Figure ||] because it simplifies the discussion later. 

What can be gathered from Figure ^| (top)? The deep valley along the diagonal is due to the fact that, provided 
that the effective SNR is fixed, we expect the BEP to be smallest when the branch power is balanced (5i = S2, 



a = 1) QDDJS02[ ]. Somewhat less expected were (a) the wide tolerance region where |5i — S2I is large (up to 
12 dB) but the BER is still small, (b) a very sharp decline in performance at the edge of the tolerance region, and 
(c) a region of high local variability in the upper part of the diagonal. The surface is truncated at 

min {b k \ = 3.75 x 1(T 6 

l<fc<M 

because smaller estimates of the (expected) BEP require an enormous computation time due to the convergence 
properties of Monte Carlo Estimation (more on this below). 



3.1 Statistically Significant Sampling Methodology 

The initial sample looks reasonable and uncovers interesting trends in system performance, but it does not contain 
enough information to make statistically significant claims. Estimating the probability that a configuration exhibits 
acceptable average performance requires several samples per point c^. The simulation is computationally expensive 
and different regions of the configuration space exhibit different variability. Therefore, we must define tight stopping 
criteria for sampling. Figure || (bottom) shows the output obtained with the following (per point c^) stopping criteria. 
The criteria are designed to achieve high estimation accuracy. 
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0.0004 0.0006 0.0008 0.0010 0.0012 0.0014 0.0016 0.0018 

BER 



Figure 5: Empirical CDF of 21 samples for a randomly chosen point vs. that of the Gaussian distribution with 
appropriate mean and variance. 

1. Sampling {x k j} stops when the relative error in the estimate b\ of the expected BEP E\b k ] is smaller than the 
relative accuracy threshold (3 = 0.1 times the current estimate b\, at a 7 = 0.9 confidence level, i.e., when 

P(\E[b k ]-b k \ <(3b k ) > 7 . 

We required n k > 2 samples to obtain an estimate a\ of the BEP variance o\. Notice that the target is the 
relative error, not the absolute error, because the range of {bk}fc=i in the configuration space spans four orders 
of magnitude. Therefore, absolute error measures are misleading. 

2. Sampling {x k j} also stops when we can say, with confidence 7 = 0.9, that the expected BEP E[b k ] is below 
the sampling threshold t = 10~ 4 , i.e., when 

P(E[b k ] < t) > 7 . 

This work considers voice quality applications, so the exact value of the expected BEP is irrelevant as long 
as it is smaller than the performance threshold T = 10~ 3 . The sampling threshold t was set to an order of 
magnitude below the performance threshold T to avoid large approximation error of a fitted surface near T. 

3. Finally, sampling {x k j} stops when more than 50 samples of 10000 frames each are required to satisfy either 
of the previous rules. This rule fired in 5% of the cases, all at the boundary of the tolerance region and most 
in mid diagonal. 

Altogether, 5154 samples were collected for an average of 6.3 samples per point. Needless to say, the computational 
expense of such sampling remains too high for practical applications. While a large number of samples is typically 
desirable (for validation purposes), we will show that our data mining framework makes very effective use of data 
and thus requires fewer samples in practice. Let us now look at the data in more detail. 

3.2 Results of Statistically Significant Sampling 

It is also likely that the samples output by the WCDMA simulation are approximately Gaussian distributed. In- 
tuitively, we are simulating a large number of information bits (800000) per BEP estimate x k j, so the Lindeberg 
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condition for the Central Limit Theorem should hold. Figure || shows empirical evidence that this is the case. We 
have arbitrarily chosen one point among those with 20-30 sample values {xkj}™=i and plotted the empirical CDF of 

this sample against that of the Gaussian distribution with the mean equal to sample mean bk and the variance equal 
to sample variance a\. The curves are close to each other and the Shapiro-Wilk test yields W = 0.98 (0 < W < 1) 
and p-value of 0.88. Other points also demonstrate similar curves and high values of W, but p-values vary signifi- 
cantly. This dataset contains sufficient samples to estimate {-E^fc]}^^ with high relative accuracy, but 6.3 samples 
per point are insufficient to formally justify a Gaussian assumption. 

It is also instructive to see some measure of how the sample variance is distributed across the configuration 
space. Figures |6] and ^ show sample sizes and sample standard deviation-to-mean ratios for the samples in Figure || 
(recall that we prefer relative measures because the range 

of {h}k =1 is large). Both fi gures indicate high variance 
around the boundary of the tolerance region. This is not surprising because the edges of the tolerance region are 
relatively steep. Figure [7| also shows relatively high variance at some points inside the tolerance region. This 
is because the simulation achieved the sampling threshold t = 10 -4 and stopped before it achieved the relative 
accuracy threshold f3 = 0.1. Knowing this, one would expect a larger relative variance in the tolerance region. Let 
us examine why this is not the case. 

We treat the BEP as a continuous Gaussian random variable bk, but all sample values {xkj}™=i are discrete — they 
are ratios of two integers, the number of errors and the number of bits simulated. The simulation may not detect any 
bit errors when the expected BEP E[bk] is relatively small (e.g., one error in the number of bits simulated). Since no 
channel is perfect, zero is too optimistic an estimate for the expected BEP. Instead, we conservatively assume that at 
least three bit errors have been detected. This is why the smallest estimate b\ of E[bk] is 3/800000 = 3.75 x 10~ 6 . 
However, using any constant cutoff prevents us from estimating the variance a\. We would need to simulate a 
large number of frames to estimate a\ when the expected BEP is small. Instead, we can empirically show that 
the probability that the expected BEP is smaller than the performance threshold T = 10~ 3 is close to one. Let 
b\ = 3.75 x 10~ 6 be the sample mean, rij, = 2 be the sample size, and a\ be the BEP variance at point where two 
independent simulations detected three or fewer bit errors each. Sampling {xkj} will stop because sample variance 
is zero, so the first stopping rule applies. 

We need to show that sampling can indeed stop, i.e., that the probability that the expected BEP is below the 
performance threshold T is 

P(E[b h ] < T) « F nfc _x ( T ~ h * ] > 0.995. 

This statement can only be false when (T — bk)^/nk/(Jk < 64, or &k > 2.2 x 10 -5 , almost an order of magnitude 
bigger than the conservative estimate b\ of the expected BEP E[bj,]. This is unlikely because Figure [7] (bottom) 
shows that the sample standard deviation rarely exceeds the sample mean even by half an order of magnitude. In 
other words, we do not have accurate estimates for variance a\ in the tolerance region. However, we can still 
reasonably conclude that configurations exhibit acceptable performance in this region. 

4 The Third Level of Aggregation: Regions 

Consider a set of buckets {Ck\¥ =1 with corresponding random variables {Bf.}^ =1 . Given a number of sample 
values, the framework developed in Section || allows us to estimate the probabilities {P(E[B},] < T)}jf =1 that 
buckets {Ck}jf =1 exhibit acceptable average performance. (All arguments about buckets equally apply to points 
because a point is a special case of a bucket.) This section is concerned with finding an optimal subset of random 
variables from among {Bk}kLi- This optimal subset corresponds to an optimal region of a 2D bucket space. We 
would like to find a sufficiently large admissible region R m such that we are sufficiently confident that buckets in R m 
exhibit acceptable average performance. 
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Figure 7: Sample standard deviation-to-mean ratios for Figure |2| (bottom). The top part shows the perspective plot 
and the bottom part shows the scatter plot. 
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X, Y parameters that partition the point space into buckets 

Mx , My X and Y dimensions of the bucket space 
M = Mx x My number of buckets in the bucket space 

Dx , Dy domains of X and Y 

77(771) number of buckets in region R m 

C K t m j\ i-th bucket in region R m , 1 < i < 77(771) 

n K(m,i) number of samples in bucket C K ( m ^ 

x K ( m ,i) , y K (m,i) X and Y values for bucket C K[m ^ 

Table 2: Summary of region notation. Also see Table |T} 



There are many ways to define admissibility and we are interested in adopting a definition that is both meaningful 
in the wireless domain and permits effective data mining algorithms. Among a space of such admissible regions, we 
can define different optimality criteria and data mining then reduces to searching within this space. In this paper, a 
region R m is admissible when it has a particular type of shape. We will explore three different criteria for the mining 
of optimal regions; the algorithms and these criteria are based on the work of Fukuda et al. [ YFM97 ] and have been 
adapted to the problem of mining simulation data in this paper. 

Additional notation relating buckets to regions is introduced in Table ||. Let X and Y be two discrete parameters 
to the temporal (e.g., WCDMA) simulations such that X and Y partition the point space into disjoint buckets 
{CfclfcLi- More precisely, let X, Y have ordinal domains D X ,D Y , let \D X \ = M x , \Dy \ = My, \D X \\D Y \ = M, 
and assume that the map p : Dx x Dy —* {Ck}^ =1 is bijective. In other words, X and Y define a discrete 2D space 
of buckets. Since the domains of X and Y are ordinal, this space is easily visualized as a 2D color map or a 3D 
perspective plot. 

For example, the average SNRs Si and S2 in the previous section partition the space of configurations into 
buckets. Both Si and S 2 vary from 3 to 42 in steps of 1 (in dB), so M x = My = 40 and M = 40 X 40 = 1600 
(recall, from Section |3[ that only 820 of these points were simulated and the remaining ones were symmetrically 
reflected). Furthermore, the domains of Si and S2 are ordinal because the values of Si and S2 are directly related to 
the powers of the transmitter antennas. In this case, the buckets are simply the points in the space of configurations. 
In general, buckets can be convex combinations of points, as detailed in Section ||. Recall that we defined the color 
of a bucket as the probability that the bucket exhibits acceptable average performance. Figure [| shows these 'colors' 
as a perspective plot for the STTD example. 



4.1 Region Shape 

Consider regions (subsets) of buckets in the bucket space. If the shape of these regions is unconstrained, there 
are 2 M possible regions {i?™}^!^. Let region R rn , 1 < m < 2 M , consist of buckets {C K ( m i)}$=i , where 
77(777.), 1 < m < 2 M , is a mapping from region number m to the number of buckets in this region, and /t(m, i), 
l<m<2 M ,l<i< r)(m), is a mapping from region number m and bucket number % within region R rn to bucket 
number k, 1 < k < M, that we use to subscript buckets {Ck}^ =l . The exact definitions of 77(771) and K(m, i) are 
not important as long as they generate all possible regions (subsets) {R m } : £ 1 =1 . 

The shape of admissible regions should be constrained because unconstrained regions are hard to interpret and 
tend to overfit the training data. Besides, the problem of selecting an optimal unconstrained region is computationally 
intractable — all 2 M possible regions must be considered, where M = 1600 in the STTD example. The region shape 
can be constrained in a number of different ways (rectangular, x-monotone, etc.). Our restrictions on region shape 
are discussed next. 
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Figure 8: Probabilities {P(E[b},] < T)}%L 1 that configurations {c^\^ =1 exhibit acceptable performance with respect 
to the performance threshold T = 10~ 3 (voice quality system). This perspective plot corresponds to the STTD 
dataset in Figure ^(bottom). The axes S\ and 52 are rotated 180 degrees counter-clockwise to provide a better view 
of the surface. 




Figure 9: Some types of admissible (connected rectilinear) regions. When we look at an admissible region from left 
to right, its upper boundary must first increase and then decrease monotonically, and its lower boundary must first 
decrease and then increase monotonically. 
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Without loss of generality, assume that Dx = {1,2,... , Mj} and Dy = {1,2,... , My}. Intuitively, re- 
gion R m is rectilinear when its intersection with any horizontal or vertical line is connected. More formally, re- 
gion R m is rectilinear if and only if whenever buckets C K ( m ,j) at (^(m.i) > 2/«(m,i) ) and C s ( mj ) at (i^), 
are both in R m , then (a) /)(r, s) = C K i m ,i) and /?(r, t) = C K r m ,j) imply buckets p(r, u) are also in R m for all 
it G [s,t], and (b) p(r, t) = C K r m {\ and p(s,t) = C R i m ,j) imply buckets p(u,t) are also in i? m for all u G [r, s]. 
Here [a, 6] means all integers between the integers a, b, inclusive. We use Manhattan geometry to define connect- 
edness. Region R m is connected if and only if for every pair of buckets C K t m> i) an d C re ( m j) m Rm there exists a 
sequence of buckets 

Cre(m,i) C/t(m,Zi)i Cre(rre,i2)' ' ' ' ' K(m,l n ) ^K(m.j) 

in i2 m such that for every 1 < k < n 

II P _1 (CK(m,/ fc )) - P~ (C/t(m,J fc+ i)) Il x = 1- 

Furthermore, we say that region R m is admissible if it is both rectilinear and connected. 

This definition of admissibility can be viewed as a relaxed definition of convexity. Geometrically, it is easy to 
see that region R m is admissible if and only if, when we look at R m from left to right, its upper boundary first 
increases and then decreases monotonically (a pseudoconcave function), and its lower boundary first decreases and 
then increases monotonically (a pseudoconvex function). In other words, the region boundary need not be strictly 
convex or strictly concave, but it must be pseudoconvex or pseudoconcave. Admissible regions are informally 
summarized in Figure ^. All admissible regions are composed of regions of four primitive types: W (region gets 
wider from left to right), N (region gets narrower), U (region slants up), and D (region slants down). Twelve 
combinations of these types yield all types of admissible regions: W, WU, WUN, WD, WDN, WN, UN, DN, U, D, 
N, and the empty region. 

Our choice of connected rectilinear regions is due to primarily heuristic considerations. These considerations 
are commonly applicable, but must be re-evaluated for each study. Both the connectedness and the rectilinearity 
restrictions can be justified for the STTD example (see next section). In general, it is easy to justify connectedness, 
but hard to justify rectilinearity. We advocate the use of connected rectilinear regions primarily because this shape 
is resistant to noise in the sample, not because we can analytically show that the region boundary is rectilinear. 
In data mining, the choice of region shape is most commonly dictated by the desired tradeoff between bias and 



variance [ ]HTF01| |. Regions with flexible shape exhibit small bias (they can fit any data) but high variance (they can 
be overly sensitive to a particular dataset). Regions with rigid shape exhibit high bias but small variance. Connected 
rectilinear regions provide a reasonable tradeoff between bias and variance for many applications. 

4.2 Evaluating Regions 

Another prerequisite to finding regions with the desired properties is a definition of region 'goodness'. Let us map 
bucket confidence P{E[B K ^ m ^] < T) to a discrete range [0 . . . 1000] and define the hit of bucket C K ( m>i } as 

h K{m>i) = [W00P(E[B K(m ^] <T) + 0.5J 

([X\ denotes the largest integer that does not exceed X), the support of bucket C K r m A as 

S«(m,i) = 1000 

(this constant was chosen to make the discretization error reasonably small), the hit of region R m as 

•q(m) 
i=l 
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and the support of region R m as 

r)(ra) 

S m =Y^ s n{m,i) = 1000r/(m). 

i=l 

The key to efficient computation of optimized-confidence and optimized-support admissible regions is the definition 
of region confidence as 

where H m is the hit and S m is the support of region R m . Let us explore the implications of these definitions in more 
detail. 



4.2.1 Model-Based and Model-Free Analyses 

Suppose, n K r m A = 6 samples have been collected for bucket C K i m j\ that consists of a single point. Let the sample 
mean be B R ( m n = 5 x 10~ 4 and the sample standard deviation be fl K r mji \ = 8.87 x 10~ 4 . Furthermore, suppose 
that five of these samples have the BER below 10~ 3 and one has the BER above 1CT 3 . Then, 

A purely model-free approach would interpret the above simulation results as 'bucket C K c m ,i) W11 l exhibit acceptable 
average performance in 5 out of 6 cases.' A strongly model-based approach would interpret the simulation results 
as 'we are 88.7% confident that bucket C K r m j\ exhibits acceptable average performance.' Our interpretation lies 
between the model-based approach and a model-free approach and posits that 'bucket C K i m>i \ will exhibit acceptable 
average performance in 887 out of 1000 cases.' These interpretations provide confidence estimates under different 
simplifying assumptions. 

The model-free interpretation does not take either sample variance or sample distribution into account. This 
interpretation is only reliable for a sufficiently large number of samples, which is a luxury in our application. Our 
middle-ground interpretation explicitly accounts for sample variance and sample distribution. When sample size is 
small, our interpretation provides a statistically valid estimate of confidence that the bucket exhibits acceptable aver- 
age performance. For a single bucket, this interpretation is as good as a strongly model-based interpretation, modulo 
a reasonably small discretization error. However, our interpretation diverges from the model-based interpretation at 
the region level. 

A strongly model-based analysis procedure would define a region random variable 

n(m) 

W m f"f K ' ^ ' 

where {B K ^ m a}^^^ are bucket random variables, {w K ( mj i)Yl=i are a priori (positive) constant weights, and 

W m = ^ W K(m,i) 
i=l 

is a normalization factor that maps these weights to probabilities of bucket occurrence in the region. A procedure 
similar to that in Section || would then be used to estimate P(E[Q m ] < T) for a threshold T. The result of this 
calculation can be interpreted as the probability that region R m exhibits acceptable average performance, condi- 
tional on the temporal simulation assumptions, the bucketing prior probabilities, and the region prior probabilities. 
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However, as we shall see later, this definition of region confidence violates a property that permits an efficient data 
mining algorithm. 

We think of region confidence in terms of average bucket confidence over the whole region, namely, 

T)(m) 

O m « — - V P(E[B K{m:i) ] < T). 
rj[m) ~ 

(If region size rj{m) is large enough, we can reasonably expect the discretization errors to cancel each other.) This 
interpretation of m does not correspond to the strongly model-based probability that region R m exhibits acceptable 
average performance. Instead, we define a region random variable P m as the probability that any bucket C K i mti \ 
in region R m exhibits acceptable average performance. Then, we estimate the expected value E[P m ] across the 
region R m by the sample mean P m Q m of estimates of bucket confidences {P(E[B K ( m a) < T)Y£l^ . 

How do these two definitions relate to each other? It is easy to show that they are equivalent only under very 
restrictive assumptions. Basically, we are assuming that the buckets are mutually independent, that population 
variance is small, and that the region is consistent, i.e., 'good' and 'bad' buckets are never mixed in the same region. 
Let bucket random variables {B K ^ m ^}^^ be mutually independent, let the estimates {S^ mi s}^ij of bucket 

variances be approximately equal to zero, and let the estimates {B^ m o}^ 1 of bucket expected BEPs be either all 

greater than the performance threshold T or all smaller than the performance threshold T (i.e., all {T — B K r m ^}^^ 
have the same sign). Then, bucket confidences 



P(E[B K{m>i) ] < T) « F nR 



rp TD 

1 - D K(m,i) 



(m,i) 



1 < i < 7](m), will be either all approximately equal to zero (B K ( m ,i) > T), or all approximately equal to one 
(B 

K(m,i) ^ T)- Therefore, region confidence @ m will be approximately equal to zero or one. Likewise, the strongly 
model-based region confidence 



P(E[Q m ] < T) w F v{m y 



Qr 



*m/y/ri(rn)J 

will be approximately equal to zero or one because the estimate ^ 2 m of region variance is (see Section ||) 



n(m) 

m i=l 



The sign of T — Q m determines whether P(E[Q m ] < T) is approximately equal to zero or one. After a minor 
rearrangement of terms, 

T - Q m = rr— W K (m,i){T ~ B n{m,i))- 



W„. , 

1=1 



We assumed that {T— B K i m ,i)}^i have the same sign, so we have shown that P(E[Q m ] < T) « Q m . The equality 

]2 1 

J K,(m,i) J 



is asymptotically exact as all variance estimates {S^ i s}^ 1 approach zero. This argument applies regardless of 



the distributions of {B K ^ m i ^ j } 1 f^ > , as long as these random variables are mutually independent. 
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4.3 Optimized Regions 

We now pursue the definition of optimized regions. Given a slope t, < r < 1, define the gain of region R m , 

1 < m < 2 M , as 

G(R m , r) = H rn — TS rn , 

where H rn is the region hit and S m is the region support. Let an optimized- gain admissible region R T with respect 
to slope t, < r < 1, be an admissible region with the maximum gain G(R t ,t) over all admissible regions 
(this region need not be unique). Optimized-gain admissible regions are easy to define, compute, and analyze, but 
hard to interpret. Common practice is to define optimized-confidence and optimized-support admissible regions. 
Admissible region R* is an optimized-confidence admissible region with respect to a given support threshold 100O?, 
< f] < M, if R* has the maximum confidence 0* = among all admissible regions with support of 

at least 1000?7. Likewise, admissible region f? is an optimized-support admissible region with respect to a given 
confidence threshold 9, < 9 < 1, if R<> has the maximum support 5<> = 1000r/(o) among all admissible regions 
with confidence of at least 9. In other words, we can either fix the region confidence 6 and find the largest region R 
with confidence of at least 9, or we can fix the minimum region size (support) lOOOr/ and find the most confident 
region R* with support of at least 1000??. 

Observe that r in the definition of an optimized-gain admissible region is the relative importance of support vs. 
that of confidence. We can find a small region with high confidence or a large region with small confidence, but both 
objectives cannot be maximized simultaneously. Increasing r will increase the confidence of the optimized-gain ad- 
missible region, but decrease its support. Likewise, decreasing r will decrease the confidence of the optimized-gain 
admissible region, but increase its support. Therefore, we can find approximate optimized-confidence and optimized- 
support admissible regions by a binary search for the value of r where the respective threshold is barely satisfied. 



The search can stop at a given level of precision At, where the lower bound on At can be found in [ ]YFM97| ] 
(they show that the number of steps in this search is logarithmic in the support 1000M of the bucket space). This 
algorithm is approximate because an optimized-confidence (resp. optimized-support) admissible region need not 



be an optimized-gain admissible region for any value of r. Yoda et al. [ |YFM97| ] argue that this approximation is 
reasonable for large datasets. 

Let us revisit the definition of region 'goodness'. Geometrically, the buckets with the same value of X are 
the columns and the buckets with the same value of Y are the rows. An optimized-gain admissible region can be 
computed in O(MxMy) time by a set of rules of the following form. Recall that a region of type W gets wider 
from left to right (see Figure |J). Let Rw(m, [s, t]) be the region of type W with maximum gain f\y(m, [s, t]) over 
all admissible regions of type W that end in column m and span rows s through t in this column. Then, either (a) m 
is the first column of R\y(Tn, [s, t]), or (b) Rw(m, [s, tj) includes the region Ry/{m — 1, [s', t']) with the maximum 
gain fw{ m , W, t']) over all admissible regions that end in column m — 1 and span rows s' > s through t' < t in 



this column. JYFM97[ ] keeps the regions with maximum gain for every region type and every triple (m, [s, t]) in a 
dynamic programming table. These locally maximal regions then grow according to a set of rules that compute an 
optimized-gain admissible region. This efficient greedy algorithm for computing optimized-gain admissible regions 
depends on the property of the gain function that we refer to as monotonicity. Let < r < 1 be a slope and R m i 
and R m ii be two admissible regions with gains G{R m i , r) > G(R m " , r). The gain function G(R m , r) is monotonic 
if for any region R^ disjoint with both R m / and R m ii 

G{R m , UR k ,T)> G(R m i> UR k ,T), 

where the union of regions is defined in the obvious way. It is easy to see that our gain function 

rj(m) 

G(R m ,T ) = H m - TS m = lW00P{E[B K(m<i) } <T) + 0.5J - IOOOt^itl) 

i=i 
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is monotonic because it is additive. However, a strongly model-based gain function 

G^ M \R m ,T) = P{E[Q m ] < T) — r V (m)/M 

is not monotonic even if we assume independence of bucket random variables {B K r m ,i)}i^ that make up Q m . To 
the best of our knowledge, only monotonic gain functions are known to result in practical algorithms for computing 
optimized-gain admissible regions. 

What happens when no estimates of mean and/or variance are available for some bucket C\P. The answer to this 
question depends on problem-specific considerations. As was demonstrated in Section ^ it is sometimes possible to 
provide conservative estimates for these values. For example, we have empirically shown that the expected BEPs of 
some configurations {c^} are smaller than T = 10 -3 with confidence P(E[bf c ] < T) > 0.995. Likewise, we know 
that as the effective SNR approaches negative infinity (in dB), the BEP approaches 0.5, which is the probability of 
correctly guessing the value of a random bit when the transmitter is turned off. Thus, we can let P(E[bk] < T) = 
for points with sufficiently small effective SNRs and a reasonable performance threshold T. If no such estimates are 
available, we can simply omit the missing buckets from the probability computation. This must be done with care 
because such buckets will contribute nothing to the confidence of the region. This fact can be used to reduce the 
computational expense of sampling. 

This section has highlighted the sometimes contradictory objectives that aggregation must satisfy: permit valid 
statistical interpretations and afford structure that can be exploited by data mining algorithms. Our approach has 
been a judicious mix of concepts from both statistics and data mining. We showed that our formulation of the data 
mining problem lies between the completely model-free approach and the strongly model-based approach. The next 
section applies the data mining methodology described here to the example in Section [| 

5 Optimized-Support Regions for the STTD Example 

This section continues the example in Section [3j. First, we show that optimized-gain regions are both rectilinear 
and connected for this example. It immediately follows that optimized-support and optimized-confidence regions 
are also admissible. An optimized-support admissible region is presented next. We show that the elaborate region 
mining setup leads to simple engineering interpretations. Finally, we look at the performance of data mining when 
the number of samples is small. Three-fold cross-validation shows that data mining performs well under these 
circumstances. 

5.1 Justification of Data Mining for the STTD Example 

Let the average SNRs S\ = X and S2 = Y partition the space of configurations in Figure ^ into disjoint points 
(buckets) {ck}^L v 1 < M < 1600. We now give an intuitive argument to justify the suitability of the data mining 
algorithm for the STTD study. Without loss of generality, consider only the points with X < Y, i.e., S\ < S%. It is 
easy to extend all arguments to X > Y, but this adds little to the discussion. 

Let c\ at (xi,yi) and C2 at (x2,yi), x\ < X2 < yi, be two points in an optimized-gain region (of arbitrary 
shape) for some slope < r < 1 (see Figure [R]). This means that the confidences of these points are one, and thus 
the expected BEPs of these points are smaller than the performance threshold T. When x±,X2 < y\ and y\ is fixed, 
the BEP is a monotonically decreasing function of x — increasing x decreases the power imbalance and increases the 
effective SNR, so the BEP must decrease. Therefore, the expected BEP of any point c u at (x u , y\), x\ < x u < X2, 
is below the performance threshold T. Thus, the confidences of points {c u } are one and these points must also be in 
the optimized-gain region R T . Three more symmetric arguments of this kind show that optimized-gain regions are 
rectilinear. 
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Figure 10: Points for arguments about region shape (see text). 



Likewise, let c\ at (xi, yi) and C2 at (x2, V2),xi < £2 < V\ < 2/2, be two points in an optimized-gain rectilinear 
region (refer to Figure [Kj. Since c\ is in the optimized-gain region and x\ < x-i, the point at (x2, yi) is also in this 
region because it has a smaller BEP than c\. Since the optimized-gain region is rectilinear, there is a horizontal path 
from (x\, y±) to (X2, yi) and a vertical path from (x2, yi) to (X2, 2/2)- Thus, there is a Manhattan path from (x\, y\) 
to (x2, 2/2)- Arguments of this kind show that optimized-gain rectilinear regions must be connected as long as they 
are 'wide enough'. 

To summarize, we have shown that optimized-gain (and thus optimized-support and optimized-confidence) re- 
gions are admissible. The data mining algorithm described in Section ||, which results in optimal admissible regions, 
is thus appropriate for the STTD example. We now show and interpret data mining results. 

5.2 Optimized-Support Admissible Regions 

Figure [TT] shows an optimized-support admissible region for the confidence threshold 9 = 0.99. Intuitively, this 
is the largest admissible region where we can claim, with confidence of at least 0.99, that configurations exhibit 
acceptable performance. This claim is conditional on temporal simulation assumptions and on mutual independence 
of configurations in the region. The shape of this region confirms that, under a fixed effective SNR, the BEP is 
minimal when the average SNRs of the two branches are equal. The width of this region shows the largest acceptable 
power imbalance. For this example, the system tolerates power imbalance of up to 12 dB. However, the width of the 
optimized region is not uniform. The region is narrower for small effective SNRs and wider for large effective SNRs. 
This means that configurations with low effective SNRs are more sensitive to power imbalance than configurations 
with high effective SNRs. None of these observations are news to an informed reader. The contribution of data 
mining in this context is not qualitative discoveries; it is statistically significant quantitative results. 

Let us see how the data mining algorithm performs when data is scarce. The initial sample of the configuration 
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Figure 11: Optimized-support admissible region for data in Figure ||| (bottom) with the confidence threshold 9 = 0.99 
and the performance threshold T = 10~ 3 . 

space in Figure || (top) contains one sample value per bucket. The statistically significant sample in Figure |2| (bottom) 
contains at least two additional sample values per bucket (recall that we required at least two sample values to 
estimate bucket variance af). Therefore, three-fold cross-validation is the most elaborate cross-validation procedure 



that this dataset affords. The regions in the top left, bottom left, and top right of Figure 12 have been computed 
for sample values in Figure ^| (top) and the first two sample values per bucket in Figure || (bottom). Each of these 
regions has been computed with two out of the three sample values per bucket. The region in the lower right has 
been computed with all data in Figure § (bottom). All four regions are optimized-support admissible regions with 
the confidence threshold = 0.95. The regions are overlaid on top of the color-coded bucket confidence values. Red 
(dark) corresponds to low confidence and white (light) corresponds to high confidence that configuration exhibits 
acceptable average performance w.r.t. the voice quality threshold T = 10~ 3 . 



The regions in Figure |12| are identical except for the lower left corner. This is not surprising because this part of 
the configuration space exhibits high relative variance. Also, the data is symmetric but the regions are asymmetric in 
the lower-left corner. Recall that optimized-gain admissible regions, and thus optimized-support admissible regions, 
are not unique. The ties in region gains are broken arbitrarily. Therefore, region asymmetry is an additional indicator 
of region instability. 



Figure |12| also shows that additional data improves image contrast but does not significantly affect region shape. 
Collecting additional sample values separates the points into ones with low confidence and ones with high confi- 
dence. A curious side effect occurs when the difference in confidence estimates of low-confidence points falls below 
the discretization error (1/1000). In this case, the 'confidence slack' 1 — 9 is allocated to arbitrary points with low 
confidence. One way to correct this situation is to raise the confidence threshold 9 — after all, more accurate data 
should afford stronger claims. Another alternative is to lower the discretization threshold. In general, optimized 
regions work best when the data is noisy. A contour plot will suffice when the data is highly accurate. 
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Figure 12: Cross-validation of optimized-support admissible regions with the confidence threshold 9 = 0.95. The 
regions in top left, bottom left, and top right have been computed with = 2 independent samples per bucket. 
There are 758 ± 2 buckets (47% of all data) per such region. The region in the bottom right has been computed 
from the statistically significant data in Figure |2| (bottom). It consists of 766 buckets (48% of all data). Red (dark) 
corresponds to low bucket confidence and white (light) corresponds to high bucket confidence w.r.t. the voice quality 
threshold T = 1(T 3 . 
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It can also be seen that the high contrast created by the sharp edge of the tolerance region is advantageous to 
data mining. The region is stable where the contrast is high. When the image is blurred, data mining tries to avoid 
the questionable boundary points. 

To summarize, this section has demonstrated that optimized-gain regions are rectilinear and connected for a 
non-trivial space of wireless system configurations. We have also shown that optimized-support admissible regions 
are easy to interpret. Finally, we have shown that data mining works well when sample sizes are small. 



6 Discussion and Future Work 



We have demonstrated a hierarchical formulation of data mining suitable for assessing performance of wireless sys- 
tem configurations. WCDMA simulation results are systematically aggregated and redescribed, leading to intuitive 
regions that allow the engineer to evaluate wireless system configuration parameters. We have shown that the as- 
sumptions about region shape and properties made by data mining algorithms can be valid in the wireless design 
context; the patterns mined hence lead to explainable and statistically valid design conclusions. As a methodol- 
ogy, data mining is thus shown to be extremely powerful when coupled with statistically meaningful performance 
evaluation. 

This work is the first (known to the authors) application of data mining methodology to solve problems in 
wireless system design. Therefore, a large number of extensions are possible and called for. We outline possible 
extensions at the three levels of aggregation: points, buckets, and regions. 

At the point level, it may be advantageous to model temporal simulations more precisely. This paper assumes a 
'large enough' number of frames per simulation and works with the distribution of estimated BEPs. We have shown 
reasonable analytical and empirical evidence that this distribution is Gaussian. The advantage of this problem formu- 
lation is the independence of spatial aggregation from the assumptions of temporal simulation. This helps introduce 
wireless engineers to the methodology of data mining for studying design problems. However, a stronger model of 
temporal simulation (e.g., Markov chains in [ Wan95T j) may yield appreciable gains in software performance. This 
direction is worth pursuing because few research groups have access to parallel computing facilities of the scale 
used in this work. For instance, the initial sample of the configuration space in Figure ^| (top) would take one year of 
computation time on a modern workstation. The study presented in this paper would clearly be impossible without 
significant computational power. 

Aggregation of points into buckets is the least developed part of this work. Suppose that we would like to simu- 
late the effects of interference on configuration performance. Assume that the distribution of the average strengths of 
the interfering signals is known a priori (e.g., estimated by ray tracing). We can either make this distribution known 
to the temporal simulation, or, alternatively, run several temporal simulations for different strengths of interfering 
signals. The former is more accurate and computationally more efficient, but the latter is more generic and simpler 
to implement. Bucketing of simulation results with varying simulation parameters is intended to approximate the 
performance of a single device under varying conditions. This paper does not employ such bucketing but instead 
builds all the necessary kinds of parameter variation into the temporal simulation (which can be argued to be the 
right way to do it). However, bucketing may be necessary when one has to work with a given dataset (e.g., measure- 
ments). Bucket space can be viewed as a configuration space for a more complex temporal simulation. Therefore, 
an in-depth treatment of bucketing is orthogonal to the primary topic of this paper, which is data mining. 

Significant work remains to be done at the region level as well. For instance, the assumption of small variance 
could conceivably be relaxed. One can also pursue the relatively difficult task of incorporating strongly model-based 
prior knowledge into the data mining algorithm, or the somewhat easier task of applying different kinds of region 
mining algorithms to problems in wireless system design. 

Defining additional case studies is another obvious direction for future work. We have studied a relatively small 
part of the parameter space of modern wireless systems. More studies of this type must be performed to highlight 
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the merits and the shortcomings of data mining in this domain. 

Finally, the strict staging of data collection and data mining can be relaxed. One can fruitfully interleave the 
two activities and have the results of data mining drive subsequent data collection. In data-scarce domains, it would 
be advantageous to focus the data collection effort on only those regions deemed most important to support a par- 
ticular data mining objective. Methodologies for closing-the-loop in this manner are becoming increasingly preva- 



lent [RBK02]. This will also help define alternative criteria for evaluating experiment designs and layouts. 
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